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ABSTRACT 

We present our results of pulse broadening time estimates and the study of the frequency 
scaling of this quantity for 60 pulsars based on actual multi-frequency scattering estimates. 
This research was based on our own measurements, performed on the observational data and 
the prohles from various pulsar prohle databases, as well as the scatter time measurements 
that were found in the literature. We were able to construct a database of over 60 pulsars with 
true multi-frequency a measurements, which allowed us to revise the previously proposed 
relations between the scatter time spectral slope and the dispersion measure (DM). We found 
that the deviations from theoretical predictions of the value of a appear for pulsars regardless 
of their DM, however the DM-averaged value of the scaling index is almost constant except 
for pulsars with very high DMs. Based on the obtained slopes we were also able to estimate the 
amount of scattering at the standard frequency of 1 GHz. We found that while the estimated 
standardized pulse broadening time increases with DM the relation seems to be much flatter 
than it was previously proposed, which suggests higher values of the scatter time for mid-DM 
pulsars, and lower values of expected pulse broadening for highly dispersed sources. 
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1 INTRODUCTION 

The temporal and angular broadening of the pulsar profiles is 
caused by the scattering of the radio waves emitted by a pulsar. 
This, along with other effects such as an interstellar dispersion of 
the signal, as well as the interstellar scintillations, happens in the 
ionized interstellar medium (ISM). Since the discovery of these 
phenomena various analyses, both observational and theoretical, 
tried to explain the properties of the in teraction of the pulsar ra- 
diation with ISM. In our recent paper jLewandowski et ^120131 
hereafter Paper 1) we performed a statistical analysis of the scatter¬ 
ing phenomenon, especially the sc atter time frequency depend ence, 
expanding upon an earlier work of iLdhmer et al.l ( I200ll . l2003 . LOl 
and L04 hereafter). Our analysis confirmed their findings that the 
frequency dependence of the scatter time measurements for the low 
dispersion measure (DM) pulsars seems to conform the theoretical 
predictions, while higher-DM pulsars apparently deviate from the 
theory, in such way that their spectral depende nce is much fl atter. 

The scattering phenomenon theory (see iRick^ 1 1 990l . for a 
review) predicts, that even if a pulse leaves the pulsar in a single 
instant, due to the propagation in the interstellar medium and the 
scattering (bending) of the rays the observer will measure this sig¬ 
nal over a certain interval. Since the radiation that had undergone 
scattering was travelling along longer geometrical paths it arrives 

* E-mail: boe@astro.ia.uz.zgora.pl 


delayed, causing the apparent pulse to attain the so called “scatter¬ 
ing tail” which is roughly consistent with an exponential decay of 
the signal. 

If one assumes the simplest geometrical model of the ISM, 
i.e. a configuration in which all the scattering happens within a thin 
screen of ionized medium, the shape of the tail may be predicted 
by assuming a two-dimensional gaussian brightness distribution of 
the scattering. In such case the resulting pulse broadening function 
(PBF) has a simple exponential decay form: PBF(t) ~ exp(t/Tj) 
where characteristic d ecay time i s called the scatter time (or 
pulse broadening time, IScheueilll968h . 

The brightness distribution of the scattering (and in return 
the shape of the scattering tail) can be estimated by assuming 
a model of the turbulent interstellar medium that is causing this 
phenomenon. Usuall y it is done b y assuming the homogeneous 
isotropic turbulence ( lRicketll[l977l) within a certain range of the 
fluctuation scales (the inner scale kj' and the outer scale fc“*). The 
spectrum of such turbulence can be described by: 

where q is an amplitude of a three-dimensional wave number and 
is the fluctuation strength for a given line-of-sight. If <K <K 
kj then the above equation simplifies to a power-law: 

P„M = (2) 
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The fluctuation spectral index 0 has to he lower than 4 
jRomani. Naravan & Blandfordill985) and in such case the scat¬ 
ter time frequency dependence also takes a form of a power law: 
Tj oc The fluctuation spectral index p and the scatter time fre¬ 
quency scaling index a are bound by a relation: a = ip 1(^—2). One 
has to note that the scatter time frequency dependence slope is the 
same as the decorrelation bandwidth slope Av^ (apart for the minus 
sign), a quantity from the interstellar scintillation (ISS) theory, in 
which both these values are bound by a formula: 

2:/rTrfAvrf = C|, (3) 

where the Ci is a constant of the order of unity, although 
it may depend on the t urbulence model or the ISM geometry 
jLambert & Rick^l 19991) . 

If one adopts a purely Kolmogorov spectrum of the density 
irregularities {p = 11/3) the resulting value of the scatter time 
spectral index is a = 4.4. Other turbulence models may yield dif¬ 
ferent values of a, but the lowest value for a thin screen model is 
a = 4.0 (for the so called critical tu rbulence model with p = A, 
[Romani. Naravan & Blandfordll 19861) . 

As we pointed out in Paper 1 (following other authors) the 
simplification of Eq.[T]may not be valid, since we do not know what 
the inner/outer scales are. Especially the inner scale may be impor¬ 
tant, since when the diffractive scale drops below k/* (and thus 
q> kj) the dependence becomes quadratic. As the diffractive scale 
changes with the observing frequency, then for a particular pulsar 
this may result in it having the theoretical scatter time slope at high 
frequencies, and the slope will become flatter at lower frequencies 
(i.e for Kolmogorov turbulence a would fall from 4.4 to 4.0). So far 
the sparsity of multi-frequency scatter time measurements made it 
virtually impossible to measure such effects. There were however 
attempts to detect the inner scale influence by the means of th e 
scattering tail shape analysis (see for example iRickett et alj2009h . 

One has to also note that the measured slopes of the scat¬ 
ter time frequency dependence may be affected by other effects, 
the most important being the anisotropy and inhomogeneity of the 
ISM. The anisotropy may falsify single scatter measurement, since 
it may affect the shape of the scattering tail jRickett et al.ll200^ 
and falsify the results of the shape modelling (yielding erroneous 
values of tj). Inhomogeneity on the o ther hand can cause the scat¬ 
tering parameters to vary in time (see ISrisken alj20ld and also 
our Paper 1), which in turn can make any attempts of estimating 
the scatter time spectrum impossible, if the observations at differ¬ 
ent frequencies were made at vastly different epochs. 

The multi-frequency scatter time analysis (often coupled with 
the decorrelation bandwidth estimates from interstellar scintillation 
observations) was attempted many times in the past. The most sig¬ 
nifican t examples are the works of ICordes. Weisberg & Boriakofil 
( Il985l) . who analysed the data for 76 pulsars, although for only 5 
of them the authors were able to obtain proper estimates of the 
decorrelation bandwidth frequency sc aling slope based on esti¬ 
mates from 3 or more frequencies. Later iBhat et alJ ( l2004h reported 
the analysis of the observational data for 98 pulsars, although most 
of them were measured at a single frequency and they were able to 
estimate the frequency scaling slopes only for 15 sources (based on 
data from 2 frequencies), however one has to mention that in the 
same paper they were able to gather and analyse a large database 
of 371 decorrelation bandwidth and pulse broadening estimates that 
was available at the time. Simultaneously LOl and L04 reported the 
scatter time frequency scaling analysis estimates based exclusively 
on pulsars which had the scatter time or decorrelation bandwidth 
estimates at three or more observing frequencies, and were able to 


gather 27 such pulsars - both from their own observations as well as 
from the literature. In our Paper 1 we decided to expand their work, 
adding a few more sources but not all of them had measurements 
on three frequencies. Hence we expanded our database by looking 
for any scatter time measurements available in the literature. 

In our research, following LOl and L04, we focused on pul¬ 
sars which had reliable pulse broadening time estimates on at least 
three observing frequencies, which allows for a proper frequency 
scaling analysis. To increase the number of such objects we revis¬ 
ited the available profile databases - the EPN databas^H and the 
profiles obtained by Parkes telescope which can be found in the 
ATNF Databas^ - and measured scatter times for selected pul¬ 
sars. Finally we supplemented all the above by adding the re¬ 
sults from our newest GMRT and Effelsberg observing projects 
(w hich involved sea rching for gigahertz-peaked spectra pulsars, 
see lKiiak et alJlroTTI) . This yielded a number of about 100 pulsars 
with multi-frequency scatter time measurements, of which about 
40 sources (including the ones already published in Paper 1) had 
reliable measurements at three or more frequencies. In this paper 
we present the results of our study. Adding our results to the 27 
pulsars from L04 we were aiming to revise the L04 study of the 
scatter time scaling index versus the dispersion measure for 60 ob¬ 
jects in total, which to our knowledge, is the largest sample based 
on true multi-frequency estimates. These estimates allowed us also 
to extrapolate/interpolate the scattering stren gth to the standard ized 
frequency of 1 GHz (similarly to the study of iBhat et alj2004l) and 
we analysed its DM-dependence. 


2 OBSERVATIONS, DATA GATHERING AND ANALYSIS 

To extend our the scatter time measurement database published 
in Paper 1 we searched all available literature, as well as added 
our own new measurements, both based on available pulsar profile 
databases, as well as from our own observations performed with 
the Giant Meterwave Radio Telescope (GMRT, Pune, India) at the 
frequencies of 325 MHz, 610 MHz and 1170 MHz. 

The new GMRT observations were conducted between 2010 
and 2012 as a part of our ong oing project to search for the gi¬ 
gaher tz peaked spectra pulsars jKiiak et^l201 ll : IPembska et al.l 
l2014h . We used 16 MHz bandwidth of GMRT working in phased 
array mode, 0.512 ms sampling rate and 256 spectral channels (see 
iGupta et al.ll200(]|) . The main goal of this observing project was 
to measure the flux densities of selected sources (to estimate the 
shape of their spectra), hence the data was intensity-calibrated us¬ 
ing measurements of known continuous sources (like 1822-096). 
Pulsar profiles were obtained by folding of 20-30 minute integra¬ 
tions (up to 1 hour in some cases). We estimate that the sensitivity 
of these observations were of the order of 1 mjy at the frequency 
of 325 MHz up to 0.5 mJy at both 610 MHz and 1060 MHz. Such 
sensitivity was usually more than enough to measure the pulsars 
flux densities, as well as it allowed for reliable study of the pulse 
profiles. 

Another source of pulse profile shapes to analyse was the 
archived profiles that can be found in two databases in the internet, 
namely the EPN Database and the Parkes ATNF Database. While 

* http: //WWW. jb.man.ac.uk/pulsar/Resources/epn/ 

^ The ANDSATNF database, which is maintained and distributed 
by the means of CSIRO Research Data Service is available at 
https://data.csiro.au/dap/search?tn=Astronomical and 
Space Sciences not elsewhere classified 
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the latter we used mainly for the L-band profiles (1.3 to 1.6 GHz), 
the EPN database contains data from a very wide range of frequen¬ 
cies (at least for some of the sources), however some of the very old 
data were not usable due to poor signal-to-noise ratio and/or poor 
time resolution of the profiles. In our analysis we omitted such pro¬ 
files since the resulting values of the scatter times may have been 
unreliable, hence bringing more of a confusion than a valuable ex¬ 
tension of the database. 

Due to the nature of the available profile data, both from our 
own observations and for the profiles obtained from databases, we 
could not ensure the data from various observing frequencies to be 
(quasi-)simultaneous, which would be a proper way to construct the 
pulse broadening spectra (as we suggested in Paper 1). Especially 
for the profiles taken from the databases the differences between 
the epochs of the observation at different frequencies may be of the 
order of several years. This makes the data prone to the inhomo¬ 
geneity effects. This needs to be taken into account when interpret¬ 
ing the results, and it is even more of a concern for our next source 
of data. 

We have also collected the scatter time measurements from 
several publications - see the in Table in the appendix for a 
full reference list. The specifics concerning the methods used for 
obtaining the scatter time measurements can be found in their re¬ 
spective references. In some cases (which was especially true for 
the older publications) the authors did not always quote the mea¬ 
surement uncertainties. To counter that we assumed the uncertain¬ 
ties to be at the level of 30%, which is probably an overestimation 
in most cases, but we found it to be appropriate at least for some of 
the lowest frequency measurements. Another problem with some 
of the measurements is that the authors of some of the oldest publi¬ 
cations we used did not publish the profiles upon which their mea¬ 
surements were based, so we could not asses their quality. This has 
to be taken into account, and as a result we were forced to disre¬ 
gard some of the measurements from our scatter time frequency 
dependence models. 


2.1 Data Analysis and Results 

For the profiles we had available (both from our observations as 
well as from the EPN and Parkes databases) we performed the mea¬ 
surements of the scatter time using the same method as in Paper 1, 
which is basically following the method used by LOl and L04. We 
have to note that there are other methods which have been shown 
successful for finding the scatter ti me for pulsars. On e of them is 
the “CLEAN” method developed bv iBhat et alj ( l2003h which does 
not require any assumptions o n the intrinsic pul se shape and shows 
very promissing results (see [Bhat et ai1l2004h . One h as to men¬ 
tion al so the cyclic spectral analysis method proposed bv IPemoresj 
( I 2 OI ll) , this method however requires access to the raw voltage 
data we did not have. We decided however against these methods, 
since more than often we were forced to combine our results with 
older Tj measurements taken from the literature (to estimate the 
slope of the scatter time frequency dependence). We decided that it 
will be more reliable if we used the same (or at least similar) meth¬ 
ods in our scatter time estimates. Additionally LOl demonstrated 
that at least in case of higher-DM pulsars the results from their 
method and the “CLEAN” method are usually consistent within 
uncertainties. 

The apparent pulse profile is a convolution of the intrinsic pro¬ 
file shape P‘(t) with functions representing several effects that af¬ 
fect the radiation, such as the scattering smearing function in the 
ISM s{t), the ISM dispersion smearing function d(t) and the receiv¬ 


ing backend response fu nction ijt). Thus, the observ ed profile shape 
may be obtained as jRamachandran et'^ll997h : 

pO(f) = p‘{t) * s{t) * d{t) * i(t), (4) 


where the asterisks (*) denote the convolution. 

The last two factors in this equation may become insignificant. 
The modem receiver systems usually are fast enough for the effect 
of i{t) to be negligible. For the coherent dedispersion receivers, or 
for filterbank receivers with large number of spectral channels the 
dispersion smearing function should also not pose a problem, espe¬ 
cially at higher frequencies. The observations at lower frequencies 
may suffer the effect, and hence the low-frequency scatter measure¬ 
ments may be falsified, especially if the actual scattering smearing 
is similar to the dispersion smearing within a single spectral chan¬ 
nel. The dispersion smearing may also be a factor at higher observ¬ 
ing frequencies, when using filterbanks with only a few spectral 
channels (as it was the case for our measurements with Effelsberg 
Radio Telescope; see the discussion in Paper 1), and it may also 
play a role for some of the older archived profile data. 

The term in the convolution (Eq. 0 that solely interest us is 
the s{t) function, which describes the broadening of a pulse due to 
the interstellar scattering. In the past various authors used various 
mathematical forms of this function. In our analysis we followed 
the simplest approach which was shown by L04, i.e. that a sim¬ 
ple exponential decay function seems to work the best, at least for 
higher DM pulsars. This fu nction describes scatterin g in the sin¬ 
gle thin screen model (see IWilliamsot]iri972L Il9’^ . and can be 
expressed as (L04): 


PBF,(t) = exp(-f/Trf)(/(0, 


(5) 


where U(t) is the unit step function (t/(t < 0) = 0, and U(t > 0) = 

1 ). 

To find the values of for the observed profiles we performed 
a least squares fit of a model shape to the observational data, to 
find a best fit parameters by minimizing the;^'^ function defined by 
the differences between the observed and modelled profiles. The 
fit was performed using the Origin softwar^ assuming a gaus- 
sian shape of the intrinsic pulse profile, or multiple gaussians in 
case of multi-component profiles. As it was discussed in Paper 1 
such approach may be viewed as simplistic, but we do not know 
the intrinsic pulse profiles. This is a problem even if one uses the 
non-scattered high frequency profiles as templates. Such approach 
may suffer from various sources of error, such as the unknown pro¬ 
file evolution (which is cloak ed by the scatterin g smear) and the 
radius-to-frequency mapping jKiiak & Gilll2003l T In Paper I we 
thoroughly discussed the possible sources of rj estimation errors 
when using our method. The most obvious of them is certain in¬ 
trinsic profile shapes, especially a case of asymmetric profile with 
an intrinsic scatter-like tail. This may not affect the resulting scat¬ 
ter time measurements by much when the interstellar scattering is 
large (i.e. much greater than the asymmetry of the profile), however 
when moving to higher observing frequencies the measured value 
of Tj will saturate at some lowest level (which will be the char¬ 
acteristic time scale of the profile’s asymmetry) when the “real” 
scattering becomes negligible. 


^ http://www.originlab.com/ 

^ In many cases to satisfy the quality of the fit we had to allow for the 
width of the gaussian template to change in the fitting procedure, especially 
at lower frequencies. This supposedly should take cai'e of the pulse width 
vaiiations occurring due to radius-to-frequency mapping 
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Figure 1. A plot of versus the observing frequency for pulsars from our sample. The solid line on each plot is the representation of the power-law fit to 
the scatter time measurements, while the dashed line shows the slope of the thin-screen Kolmogorov prediction with a = 4.4. {Figure continued on the next 
pages). 


One has to also remember about the very large scattering 
cases, i.e. instances where the scatter time is a significant fraction 
of the pulsar rotational period, or even greater. As we discussed 
in Paper 1 in such cases the background level of the profile may be 
falsified which may lead to underestimation of the measured scatter 
time. There is however a way to salvage at least some of the infor¬ 
mation as we did in case of PSR B1557-50 observed at 325 MHz 
(and for two other pulsars), by correcting the profile background 
level. To address this we carefully examined the profiles we were 
analysing, but in our current sample we did not find observations 
which required us to use this method. However, we suspect, that 
at least some of the published measurements we used may have 
suffered from this, especially at the lowest observing frequencies. 
This can be recognized in cases where the quoted values of the 
scatter time do not seem to correspond with the values measured 
at higher frequencies. Such measurements were usually omitted in 
further analysis. 

The full listing of the new data in our scatter time measure¬ 
ment database, including our new measurements can be found in 
Table [3 (see also the online material). This data, along with the 
measurements published in Paper 1, as well as the data from LOl, 
L04 and earlier experiments (see Paper 1) constitutes the most com¬ 
plete multi-frequency measurement database compiled so far, 
and this was the object of our further study. 


3 THE FREQUENCY EVOLUTION OF PULSE 
BROADENING 

The main goal of our analysis was to construct the scatter pulse 
broadening spectra for as many pulsars as possible, to ascertain 
the frequency dependence of this quantity and compare it to the 
theoretical predictions. Using our database, which includes all the 
published data we could find as well as our new measurements (its 
acquisition was described in the previous section) we attempted to 
construct the scatter time spectra for about 80 pulsars. For the ma¬ 
jority of them it was the first attempt to do so, for several other 
sources our current analysis added some individual measure¬ 
ments and forced us to re-evaluate the previously published spectra 
- see Paper 1 (and references therein) for a summary of scattering 
measurements prior to the results of our current analysis. 

As we mentioned in the Introduction one can expect the mea¬ 
sured scatter time values to be highly dependant on the observing 
frequency, and usually following a simple power-law. If one uses a 
simple thin screen model (or for that mater also for an uniformly 
distributed medium) with Kolmogorov type turbulence, then (ex¬ 
cluding anisotropy, inhomogeneity and inner scale effects) one can 
expect the scatter time spectral index to be o- = 4.4. For other ISM 
geometries, or when the inner scale effects start to be significant 
one can expect the spectral index to drop down to a = 4.0 or lower. 
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Figure[T] A plot of vs the observing frequency for pulsars from our sample (continued). 


even with the assumption of Kolmogorov fluctuat ion distribution in 
the ISM (see for example ICordes & LaziolEoOlh . 

To find the spectral slopes of the pulse broadening frequency 
dependence we modelled our data using a power-law. For the 
proper uncertainty analysis we performed weighted least-squares 
fits of the full power-law function: 

Trf(v) = (6) 

instead of fitting a linear relation to the logarithms of and v. To 
perform weighted fits in cases where values did not have error 
estimates - which was often the case for some of the archival scatter 
time data, as well as for some more problematic of our pulse profile 
models - we assumed a 30% uncertainty. 

Figure [T] shows the results of our modelling for 48 pulsars of 
the total number of about 80 pulsars we analysed in our current 
project. The remaining objects have either been measured at two 
frequencies only, or - due to the exclusion of some of the measured 
Tj values - were limited to either two data points or points at only 
two observing frequencies. As we mentioned in the Introduction, 
following LOl and L04 we decided to omit such results in our fur¬ 
ther analysis (beyond quoting the value of a in the table) as we 
were unable to reliably estimate the uncertainties of the spectral 
index values obtained from the fits. 

We were also forced to omit some of the individual tj mea¬ 
surements from our fits as we found them unreliable. These were 
mainly the very old lowest frequency measurements (below 100 
MHz) that did not conform higher frequency data, and in most cases 


were quoted without uncertainties in their original papers. We also 
disregarded some of the high frequency measurements (including 
some of our own) where we expected that the estimate was falsi¬ 
fied by the limitations of our method - see the previous section and 
more detailed discussion in Paper 1. Using our method for profiles 
which show very little scattering yields the upper limits for the scat¬ 
ter time rather than actual estimates, and there is no reliable way to 
include these upper limits in our fits. 

Table [T] summarizes the results of our fits for the pulsars in 
our sample, including only the sources for which the scattering was 
measured on at least three frequencies (full list of the scatter time 
measurements for all the pulsars from our database can be found 
in the online appendix). In this table, apart from the basic rele- 
vant pulsar data (that w ere taken from the ATNF Pulsar Catalogu^ 
[Manchester et al.l2005l) we quote the the values of the pulse broad¬ 
ening time spectral index a, and the inferred turbulence spectral 
index f}. The fj values were calculated using the fj = la! {a - 2) 
relation, which is valid only for a > 4, hence we do not present the 
resulting beta values when this relation is irrelevant (i.e. a within 
its uncertainties does not reach the value of 4). This does not mean 
however, that the values of a that are lower than 4.0 should be con¬ 
sidered wrong or unreliable, only that within the constraints of the 
basic scattering model there is no way to estimate jj if a < 4.0, 
but such values are not “disallowed” in general (they are expected 


^ Available at http: //www. atn£. csiro. au/people/pulsar/psrcat/ 
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FigureU A plot of versus the observing frequency for pulsars from our sample (continued). 


to appear for example for scattering geometries different from the 
thin screen model). 

As one can see from the table the obtained values of 
a range from 0.73±0.15 for PSR B1054-62 to 5.61±0.50 for 
PSR B0402+61. Especially the very low values of a may be a 
source of concern and we examined all such cases closely. The most 
common reason behind such low values were the complicated mor¬ 
phologies of the profiles (in most cases clearly multi-component or 
asymmetric ones) which may have influenced our fits to the data. 
One such case was the mentioned B1920-f21 (see also Fig. [T] with 
an asymmetric profile and actually no significant scattering visible 
in any of the analysed profiles. In such cases one could say that 
our method measured not the evolution of scattering with observ¬ 
ing frequency, but the frequency evolution of the intrinsic profile. 
Similar was the reason for excluding B1845-01. Another variation 
of this are cases where we had reliable measurements at lower ob¬ 
serving frequencies but only other measurements were at very high 
frequencies where scattering was insignificant, but our method still 
gave us some values due to profile asymmetry. One such case is def¬ 
initely PSR B1054-62 (see Fig.[T]and pulse broadening measure¬ 
ments in the online appendix) for which we had data at 410 MHz 
with a reliable measurement, but the only 2 other data points are 
from our measurements performed for the profiles taken from the 
EPN Database at 1372 MHz and 1590 MHz. These profiles appear 
single but highly asymmetric and apparently the asymmetry is in¬ 
trinsic and not related to scattering. 

As we mentioned above, in our further scatter time frequency 
dependence analysis we decided to exclude all the pulsars which 
had scattering measured at only two frequencies and additionally 
for most of them the values of a are suspiciously low. We also 
excluded also several other pulsars with low values of or, like the 
pulsars which were measured at three frequencies but two of them 
were relatively close, or pulsars with multi-component profiles. All 
the or values that we omitted are quoted in italic in Table [T] We be¬ 
lieve that the remaining 33 estimates of the scatter time frequency 
scaling index or are the most reliable estimates obtained for these 
pulsars so far. 

Our final sample of pulsars with useful scatter time spectra 


estimations includes three pulsars that we already analysed in Pa¬ 
per 1, and that is because this time new profiles for these pulsars 
and/or obtained values of t from the literature that we did not use 
previously. The latter was tru e in case of PS R B1557-50, where 
we added a measurement from lRicketd ( 1 1 977l) at 750 MHz. The ad¬ 
dition of that measurement changed the a from 3.81±0.09 that was 
quoted in Paper 1 to 4.31±0.36 in our current analysis. This shows 
how sensitive are the spectral fits even to individual measurements. 
Those two values are definitely statistically different, but one has to 
note that the added measurement is from over 35 years ago while 
the remaining ones are presumably from the last few years. In such 
cases the inhomogeneity effects (see the Introduction) that are able 
to change the pulse broadening over time may have played a signif¬ 
icant role. This case strongly shows how such effects may influence 
the measured spectral slopes. 

Another pulsar published already in Paper 1 was 
PSR B1641-45. For this one the previous spectral slope was 
3.81 and was based on two frequency observations only (hence 
no error estimate) . In our curren t analysis, where we also added 
observation from iRicked (Il977h the resulting a is 3.84±0.08, 
hence in this case the 35-year old measurement appears to be 
in a very good agreement with the newer data. Finally the last 
case of the pulsar common for Paper 1 and our current analysis is 
PSR B1758-23, the source with the highest DM value (1073.9 pc 
cm“^). As a result of our previous analysis we quoted the value of 
a = 4.92 ± 0.11, while in our current approach we added new data 
obtained for profiles from the ATNF Database, as well as added an 
archival measurement from Manchester et al. (1985). The addition 
of the new data points clearly showed that the lower-frequency 
Tj measurements show a completely different slope than the high 
frequency measurements at 2.6 and 4.8 GHz. These high frequency 
data were included in the previous fit (which consisted of 4 points 
in total), but in our current analysis we had 10 measurements 
available (see the online appendix) so we were able to exclude the 
two suspicious high frequency ones. As a result the obtained a 
value is much lower and the pulse broadening slope is even even 
flatter than the results published earlier by LOl (3.91 ± 0.4, based 
on three frequencies). 
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Table 1. The scattering spectral index {a) and the electron density spectral index (j3) for the observed pulsars. Values of a quoted in italic are considered 
doubtful (see aificle text for explanation). Values of yS indicated by an asterisk are “un-physical”, as the /3 = 2a/(a — 2) relation is valid only for a ^ 4. Table 
also lists the scattering fluctuation strength log and the estimated scattering at a standard frequency of 1 GHz. In case of both of these the values given in 
brackets are considered unreliable (calculated based upon a doubtful value of a). 


Pulsar 

In 

(deg) 

bii 

(deg) 

DM 

(pc cm“^) 

Distance 

(kpc) 

a 


logCi 

Td (ms) 
at 1 GHz 

B0320+39 

152.18 

-14.34 

26.764 

1.49 

1.77±0.2 

* 

(-1.36) 

(0.29) 

B0339-t-53 

147.02 

-1.43 

67.30 

2.48 

3.31±0.2 


-3.01 

0.0174 

B0402-t61 

144.02 

7.05 

65.303 

3.05 

5.61±0.5 

3.1±0.5 

-2.86 

0.0725 

B0450-18 

217.08 

-34.09 

39.903 

0.4 

1.08±0.2 


(-0.43) 

(0.434) 

B0531-t21 

184.56 

-5.78 

56.791 

2.00 

4.21±0.5 

3.8±0.9 

-4.12 

0.00071 

B06114-22 

188.79 

2.40 

96.91 

4.74 

1.73±0.5 


(-2.07) 

(1.08) 

B0740-28 

243.77 

-2.44 


2.00 

2.52±0.3 


(-1.86) 

(0.28) 

B0808-47 

263.30 

-7.96 

228.3 

12.71 

3.17±0.16 


-1.43 

1.94 

B08094-74 

140.00 

31.62 

6.116 

0.43 

3.90±0.85 

4.1±1.8 

-4.76 

0.00025 

B0833-45 

263.55 

-2.79 

67.99 

0.28 

4.45±0.35 

3.63±0.57 

-0.98 

0.047 

B0835-41 

260.90 

-0.34 

147.29 

1.50 

2.36±0.09 


(-1.56) 

(0.353) 

B0839-53 

270.77 

-7.14 

156.5 

7.77 

4.69±0.02 

3.5±0.02 

-3.00 

0.304 

B09194-06 

225.42 

36.39 

27.271 

1.10 

2.73±0.3 

* 

(-2.70) 

(0.0077) 

B1054-62 

290.29 

-2.97 

320.3 

2.40 

0.73±0.15 


(-0.68) 

(11.4) 

B1114-41 

284.45 

18.07 

40.53 

2.68 

3.50±0.17 

* 

-1.66 

0.210 

B1154-62 

296.71 

-0.20 

325.2 

4.00 

4.59±0.27 

3.54±0.42 

-2.30 

0.493 

B1323-58 

307.50 

3.56 

287.30 

3.00 

4.00±0.08 

4.00±0.16 

-0.48 

9.47 

B1323-62 

307.07 

0.20 

318.80 

4.00 

5.57±0.55 

3.12±0.62 

-2.19 

2.35 

B1356-60 

311.24 

1.13 

293.71 

5.00 

3.77±0.39 

4.26±0.88 

-1.68 

1.04 

B1508-(-55 

91.33 

52.29 

19.613 

2.10 

1.66±0.08 

* 

(-2.37) 

(0.079) 

B1557-50 

330.69 

1.63 

260.56 

6.90 

4.31±0.36 

3.73±0.62 

-1.78 

4.64 

B1641-45 

339.19 

-0.19 

478.8 

4.50 

3.84±0.08 

4.17±0.18 

-0.54 

11.0 

B1642-03 

14.11 

26.06 

35.727 

2.91 

2.52±0.54 

* 

(-3.48) 

(0.0076) 

J1723-3659 

350.68 

-0.41 

254.2 

4.28 

3.16±0.04 

* 

0.32 

3.59 

B1749-28 

1.54 

-0.96 

50.372 

0.20 

3.78±0.08 

4.24 ±0.18 

-1.60 

0.0028 

B1758-23 

6.84 

-0.07 

1073.9 

4.00 

3.45±0.09 

* 

2.33 

306.9 

B1818-04 

25.46 

4.73 

84.435 

0.30 

2.54±0.10 

* 

(-0.16) 

(0.488) 

B1821-19 

12.28 

-3.11 

224.648 

3.70 

3.67±0.11 

* 

(-1.49) 

0.707 

B1844-04 

28.88 

-0.94 

141.979 

3.12 

4.71±0.08 

3.48±0.12 

-1.80 

0.224 

B1845-01 

31.34 

0.04 

159.531 

4.40 

2.52±0.09 

* 

(-1.49) 

(4.61) 

B18594-03 

37.21 

-0.64 

402.080 

7.00 

4.42±0.07 

3.65±0.11 

-2.30 

1.41 

B19004-01 

35.73 

-1.96 

245.167 

3.30 

3.53±0.10 

* 

-1.32 

0.58 

B19074-02 

37.60 

-2.71 

171.734 

4.50 

3.58±0.24 

* 

-3.05 

0.049 

B19074-10 

44.83 

0.99 

149.982 

4.80 

3.61±0.03 

* 

-3.28 

0.036 

B1911-04 

31.31 

-7.12 

89.385 

3.22 

2.62±0.86 

* 

-2.50 

0.141 

B19134-10 

44.71 

-0.65 

241.693 

7.00 

3.69±0.12 

* 

-2.02 

0.861 

B19144-13 

47.58 

0.45 

237.009 

4.50 

1.43±0.22 

* 

(-2.22) 

(0.637) 

B19154-13 

48.26 

0.62 

94.538 

5.00 

2.11 ±0.05 

* 

(-2.64) 

(0.249) 

B19194-21 

55.78 

3.50 

12.455 

0.30 

3.45±0.07 

* 

-3.26 

0.00031 

B19204-21 

55.28 

2.93 

217.086 

4.00 

0.94±0.88 

* 

(-1.19) 

(7.46) 

B19294-20 

55.57 

0.64 

211.151 

5.00 

4.53±0.65 

3.58±1.0 

-2.80 

0.179 

B19334-16 

52.44 

-2.09 

158.521 

3.70 

3.35±0.15 

* 

-2.89 

0.0417 

B19464-35 

70.70 

5.05 

129.075 

7.87 

3.63±0.06 

* 

-2.24 

0.619 

B19534-50 

84.79 

11.55 

31.974 

1.80 

2.59±0.01 

* 

(-2.01) 

(0.151) 

B20024-31 

69.01 

0.02 

234.820 

8.00 

4.21±0.66 

3.81±1.19 

-3.45 

0.0689 

B20534-36 

79.13 

-5.59 

97.3140 

5.00 

3.78±0.08 

* 

-2.85 

0.0964 

B22174-47 

98.38 

-7.60 

43.519 

2.45 

3.22±0.32 

* 

-4.33 

2.67 

B23034-30 

97.72 

-26.66 

49.544 

3.92 

3.42±0.26 

* 

-3.94 

8.50 


3.1 Scatter time scaling index versus the dispersion measure 

Figurel^shows the values of the scatter time spectral index a versus 
the dispersion measure. For comparison we are showing two plots: 
the re-created plot of L04 (top) and the current plot of the spectral 
slopes for 61 pulsars (L04 had 27 sources). The plot shows only 
the values of a for pulsars that had Tj measured on at least three 
frequencies (in this regard it is different from the respective plot we 
showed in Paper 1), and were not excluded for any of the reasons 


stated in above paragraphs. To create the new version of the plot 
in addition to the data used by L04 and all the references therein 
(see the figure caption for a full reference list) we used our data 
presented in Paper 1 (these are shown as open red circles) and the 
results of our current research (full blue circles). We also added 
two a estimations from our scintillation stud ies of PSR B0329-I-54 
(a = 3.86 ± 0.24. iLewandowski et alJlMH) and PSR B0823-H26 
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(a = 3.94 ± 0.36. lDaszuta. Lewandowski & Kiiakll2013l) : these are 
shown in Figurel^hy diamonds. 

The current version of the a versus DM plot looks signifi¬ 
cantly different from the one showed hy L04, and not only because 
we increased the number of measurements (by the factor of 2.25). 
First off one has to note that some of the LOl and L04 measure¬ 
ments had disappeared, as we replaced them with our own, up¬ 
dated values for the given sources. In most cases the new values 
agree with the theoretical predictions better, which in fact counters 
the LOl argument that most of the high-DM pulsars show signif¬ 
icant deviation from theoretical predictions. That is true for LOl 
measurements (open squares) for high DM pulsars, as we already 
showed in Paper 1 at least in two cases we replaced measurements 
which were way below theoretical predictions (i.e. 4.0 < a < 4.4) - 
this can be seen as two of the open squares from around DM=600 
on the top plot being replaced by two open circles with a very close 
to 4.0 on the new version of the plot. As one can see for 7 pulsars 
from LOl which had very low values of a- there is only 5 remain¬ 
ing. While the 2 (out of 7) new “theory consistent” values of a may 
not seem as much let us point out that these were the only two pul¬ 
sars for which we had new data, for the rest we had to retain the a 
measurement from LOL 

Similarly, for the mid-DM pulsars analysed in L04 there was 
a few for which the values of a were greater than 4.4 (although 
consistent with thin screen predictions within error estimates). Our 
current analysis for some of these pulsars yielded more theory- 
compliant a values, and were replaced by full circles on the middle 
plot. 

Another interesting result is that the values of a we got from 
the scattering analysis when compared to the estimations from 
decorrelation bandwidth seem to show similar spread, when we 
limit the range to DM< 300 pc cm“^. In the middle plot of Fig|2]the 
scattering-based measurements are shown as circles and squares, 
whereas the scintillation-based estimates are represented by trian¬ 
gles, diamonds and a cross. In this DM range in cases of both types 
of measurements there is a solid number of pulsars for which a 
agrees with the theory-predicted range (4.0 < a < 4.4), and the 
majority of deviations is towards lower values of a, but usually not 
much lower than ff = 3.5. 

The general idea of LOl and L04 when interpreting the a ver¬ 
sus DM plot was, that while for low-DM pulsars the scatter time 
spectral index seems to agree with theoretical predictions, the pul¬ 
sars with DM> 300 pc cm“^ deviate from these predictions signifi¬ 
cantly. In our Paper I we agreed with this , however we postulated 
that the division is still there, however it appears at slightly lower 
DM values, around DM= 230 to 250 pc cm“^. Looking at the cur¬ 
rent plot one may comment, that both the original, as well as our 
conclusion were based on data biased by a small number of scatter 
time measurements. 

The bottom plot in Fig. shows the averaged values of 
the scatter time frequency scaling index, and the averaging was 
performed for 12 log-DM bins (in a DM range between 5 and 
1100 pc cm“^). As one can see the averaged a values for the 
most part stay very close to the theoretically allowed value of 
a- = 4.0. The only notable exceptions happen for the bin around 
DM=100 pc cm“^ (which probably should be disregarded due to 
selection effects, as there are only two pulsars in that bin, and one 
of them has a suspiciously low value of a), the last bin containing 
the highest DM pulsars (DM> 700 pc cm“^) and possibly the lesser 
extent the bin before the last. There is no proof for the global devi¬ 
ation of a for pulsars with DM> 300 pc cm“^ that L04 suggested. 
In the averaged a values significant deviations from the theoretical 





Figure 2. A plot of the spectral index of pulse broadening a ver¬ 
sus the dispersion measure. Top panel shows the data from L04 
only, the middle panel includes our measurements. Empty squares 
represent th e measurements from LOl, fi ll ed squ ares - L04, empty 
triangles - | Jphn^n, Nicash' o & K oribalskil |l99j), filled Mangles - 
ICordes W eisberg & Boriak^ jT985l)7 cross - [K uzmin et alj |20 0^. dia ¬ 
mond - iLewandow sS^eT^JbOl ihJpaszuta. Lewandow.skr&Kiialdi20l^ . 
open circles (redi - iLewaiidow^ et al.l i2013l) . full circles (blue! - this pa¬ 
per. The bottom panel shows the averaged alpha values (for 12 bins in 
logDM). 


predictions seem to start at much higher values of DM. Even for the 
DM-bin covering the 500-600 pc cm“^ range the average a value is 
statistically (within the l-cr shown on the plot) in agreement with 
the theoretical predictions. As for the last DM-bin (DM> 700) we 
have to point that the deviation may be possibly be due small num¬ 
ber statistics, since its average a value is based on only three pul¬ 
sars. 

While the averaged values seem to be largely consistent with 
the simple scattering theory, for the individual pulsars large devi¬ 
ations are much more common. Apparently these start already at 
very low values of dispersion measure, and continue over the whole 
range of available DMs. If anything it is only the range of the devia- 
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Figure 3. Pulsars with measured scattering spectral index in the galactic 
X-Y coordinates with the spiral arms positions based on NE2001 model 
iCordes & Lazidl2002t) . Symbol (and colour) coding the same as in Fig|2] 
The size of the marker indicates the deviations from theoretical predictions. 
The largest markets represent pulsars with the a values below 4.0, the small¬ 
est ones - objects with a > 4.4. 


tion that seems to be increasing for the larger DMs. With 60 pulsars 
in the plot the selection effects will still be pretty strong, but our in¬ 
terpretation is that with the increasing dispersion measure (and thus 
the distance) the probability of deviation also increases, but even 
amongst the highest-DM pulsars there are some pulsars for which 
the agreement with theory is reasonably good. 

We also plotted the pulsars from out sample in the Galactic 
X-Y coordinates in Fig.[3 On this plot we present the positions of 
the pulsars, and the size of the graph marker represents the value 
of the pulse broadening spectral index. The medium sized markers 
were used for sources for which their 4.0 < a < 4.4 (within error 
estimates), the largest markers represent objects with a definitely 
lower than 4.0, and the smallest ones - pulsars with a > 4.4 (which 
is only 3 sources, corresponding to 3 smallest size full circles). 

We can safely state that there is no correlation between the 
value of the pulse broadening spectral slope and the distance and/or 
the position of the source in our Galaxy. One thing to remember 
when interpreting this plot is that the pulsar positions, especially for 
the more distant objects were based mainly on t he galactic electron 
density model (NE2001. lTavlor & Cordesll993h . and in many cases 
the precision of such distance estimation, and hence the precision 
of the actual pulsar position in the Milky Way is relatively low. 


3.2 Search for the evidence of inner scale effects 

As in Paper 1 we also searched our data for any evidence of the in¬ 
ner scale effects influencing the measurements in our database. 
Table shows the results. In case of a few selected pulsars, for 
which we had a large number of rj measurements we split the fre¬ 
quency range into two sections. The frequency at which the divi¬ 
sion happened was different for all the sources, and decided upon 


Table 2. Pulse broadening spectral slopes calculated in different frequency 
ranges for selected pulsars from our database. a/p„, is the spectral index 
in the low-frequency range, - in the high frequency range and Vrf,v 
denotes the frequency at which the data was divided. 


Pulsar 

^low 

Vrfh(MHz) 

^high 

B0809-r74 

1.51+0.28 

44 

3.90+0.85 

B1154-62 

3.79 

685 

5.24+ 0.07 

B1557-50 

4.28+0.44 

900 

3.76+0.41 

B1758-23 

3.45+0.09 

1600 

4.76+0.62 

B1821-19 

5.2+1.9 

500 

3.86 

B1859+03 

4.42+0.07 

925 

6.0+.0.6 

B1913+10 

3.66+0.13 

925 

4.47 

B1946+35 

3.62+0.06 

430 

4.00 


the view of the data (see Fig.[T]l. In selected cases we used mea¬ 
surements that were previously excluded from the overall analysis. 

As it was mentioned earlier in the Introduction, the inner scale 
effects should manifest themselves as a flattening of the spectral 
slope at low observing frequencies. One thing to mention is, that 
the frequency at which it should happen actually depends on the 
inner scale value, and it is impossible to know it beforehand. 

Again, as it was in case of our similar analysis in Paper 1, the 
results presented in Table are highly inconclusive, mainly due to 
the quality of the data and subsequently - the reliability of the scal¬ 
ing index fits, which are often based upon just two measurements. 
For instance, in the cases of B0809-)-74, B1154-62 and B1859-1-03 
one of the slope measurements is highly unrealistic. For B1557-50 
(similarly to what we noticed in Paper 1) and B1821-19 the lower 
frequency slopes are actually greater than the high frequency ones; 
an opposite to what we would expect from the inner scale effects. 
In the remaining cases - B1758-23, B1913-)-10 and B1946-I-35 - 
we got the desired effect, but especially the high frequency slopes 
are highly uncertain - in case of the first pulsar there is a large error 
bar (fit based on just 3 widely spread data points) and in case of the 
other two objects we have error-less values based on just two fre¬ 
quencies. Also, one has to note, that for a number of sources shown 
in Table|2]the slope differences are not significant statistically (usu¬ 
ally below 3-(T level). 

The above clearly shows that to perform this kind of analy¬ 
sis one needs much better quality data, especially when it comes 
to the frequency coverage of the pulse broadening measurements. 
Also, such data should be free of the inhomogeneity and anisotropy 
elfects, preferably the measurements should be performed quasi- 
simultaneously at all the frequencies, at least in the sense of the 
inhom ogeneity scale whic h is presumed to be of the order of a few 
weeks jRickett et al.ll2009l) . 


4 THE STRENGTH OE SCATTERING 

In Table [T] one can find the estimated fluctuation strength for the 
pulsar’s line-of-sight . Following L04 (see their Eq.6) we cal¬ 
culated this value using the parameters of the fit that yielded our 
scatter time spectral slope a. Using these parameters we extrapo¬ 
lated the scatter time to a uniform frequency of 1 GHz. Next, using 
the value of Ci = 1.16 and Equation [3 we calculated the expected 
decorrelation bandwidth at this frequency, which we entered into 
L04’s Eq. 6 to obtain the value of Cl . 

In most cases the results are much higher than the canoni¬ 
cal value expected for homogeneous medium (logC^^ = -3.5, see 
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Johnston. Nicastro & Koribalskj||l9^ . and this was also reported 
earlier by other authors (see for example L04). The lowest value 
we found was for PSR B0809+74 with log = -4.72 which is 
not uncommon for nearby pulsars (and this source has the lowest 
DM in our sample). The highest value we found was log = 2.33 
for B1758-23 - the highest DM pulsar. When looking at the values 
of the fluctuation strength one has to remember that the estima¬ 
tion is very sensitive to the value of the scatter time spectral index, 
especially if one underestimates a the values of may be overes¬ 
timated. Hence in cases where we are not certain of a (see previous 
paragraphs) in Table[T]we quoted the log in brackets. 

The value of the fluctuation strength is one of the ways to 
describe the strength of scattering happening along the line-of-sight 
of an object, but by definition it is supposed to be distance inde¬ 
pendent. To estim ate the actual amo unt to scattering we decided to 
follow the idea of iBhat et~^ ( l2004h and also L04. To compare the 
measured values of pulse broadening between different sources we 
re-scaled the scattering of the pulsars from our sample to a standard 
frequency of 1 GHz, using the actual frequency scaling power-law, 
i.e. the values of a and b we got from our fits (see Eq. [6}. We also 
re-calculated the values of for the frequency of 1 GHz for the 
data we published in Paper 1, as well as for all the data from the 
references used to create Fig.|2(see the figure caption for the refer¬ 
ence list), making our own fits to the scatter time or decorrelation 
bandwidth data where necessary. For the scintillation data we inter¬ 
polated (or extrapolated) the value of the decorrelation bandwidth 
to the observing frequency of 1 GHz, and then calculated the cor¬ 
responding pulse broadening time at this frequency using formula 
from Eq.O assuming (initially) the value of Ci = 1. The result is 
shown in the top panel of Figure |4]. For the purposes of this plot 
we used only objects in which we are certain about the values of 
the scatter time scaling index a (i.e. the values of at 1 GHz that 
in Tab.[T]are not quoted in brackets). 

One thing to notice about this plot is the apparent discrep¬ 
ancy between the scattering-based values of at 1 GHz and the 
scintillation based estimates. The latter are presented as diamonds 
and both types of triangles, while the remaining markers (full and 
empty circles and squares) are from scattering measurements. The 
discrepan cy is visible especially for the estim ations based upon 
data from iCordes. Weisberg & Boriakoffj ( Il985h (full triangles). It 
seems that the bunch of scattering measurements does not follow 
the trend defined by these data, and the scattering estimates are 
slightly higher on the plot or the scintillation estimates are too low. 
Such assumption may be of course biased by a relatively poor over¬ 
lap between these two datasets, and the apparent discrepancy is 
based only on a very few estimates, however we believe it to be 
true. 


To asses the discrepancy we looked at the pulsars for 
which we were able to get scatter time at 1 GHz estimates 
from both the scattering and the scintillation observations, and 
we found 4 such sources. We had to exclude 2 of them from 
start, since we disregarded our a values for B0329+54 and 
B0823-(-26. In the remaining two we have the Vela pulsar 
(B0833-45) for which our estimate from scattering measurements 
was 0.0466 m s, while based on decorrelation ban dwidth mea¬ 
surement from ICordes, Weisberg & Boriakofil l ll985h we got tj( 1 
GHz)= 0.00804 ms which means difference by a ratio of 5.8. The 
other one is B1933-t-16 for which we obtained the standard fre¬ 
quency Tj of 0.0417 from scattering and 0.00883 from scintillation 





Figure 4. Three plots of the 1 GHz scattering versus the dispersion 
measure for two values of the Ci constant: Ci = 1 (top panel) and 
Cl =5 (middle and bottom panels). The scattering-based estimates 
are shown by empty squa res - LOl, filled squares - L04, open circles - 
iLewandowski et afl i2013h . full circles - t his paper. The scintillation-based 
estima tes are shown as: em pty triangles -IJ ohnston. Nicastr o & Koribalskil 
h998h. filled t ri angles - ICordes. Weisberg & Boriakofil 1 198 A , cross 
- [ Kuzmin et alj i20Q2h. and diamo nds - ILewandowski r izoTTI) : 
iDaszuta. Lewandowski & Kiiak]i2013h . The bottom panel also includes the 
lines corresponding to our fits (the blue solid and dashed lines) as well as 
the previously proposed ver sus DM relations: the red dash-dotted line 
for iRamachandran e^^j ^1997 ), the dotted red line for L04 and the three- 
dot-dashed line for iBhat et 20041) . 
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data (again based on ICordes. Weisberg & Boriakoffll98^ - a raticH 
of 4.7. 

The difference between these values can be explained assum¬ 
ing that the Ci constant is in reality not close to the unity, hut rather 
close to 5. We explored this hypothesis and used the value of Ci =5 
to create the plot shown in the middle panel of Figure|4] One has to 
remember however that this value is based on just two pulsars, and 
the scatter time and scintillation measurements were not performed 
simultaneously. Also, our way of estimating C| by comparing the 
decorrelation bandwidth estimate and the scatter time estimate af¬ 
ter extrapolating them both to the frequency of 1 GHz may be a 
source of additional error due to the uncertainties in the scaling in¬ 
dices. Since the scaling indices for both of these quantities are not 
the same - the resulting Ci value will depend on the choice of the 
standard frequency. It would be probably to far to suggest at this 
point that the value of Ci may be frequency dependant (that would 
require us to re-evaluate the whole scattering and scintillation the¬ 
ory), but in our approach that would be an expected outcome. 

Finally, the pulsars we used are in some regards fairly non¬ 
standard: the Vela pulsar profile shape makes it extremely hard to 
obtain a reliable scatter time measurement (tbe intrinsic pulsar pro¬ 
file apparently exhibits a scatter-like tail), whereas PSR B1933-I-16 
shows an extremely low scaling index {a = 3.35) and the peculiar¬ 
ity of its line-of-sight was discussed by many authors before (see 
L04 for example). 

The value of Ci we suggest h ere is also muc h larger than unity 
which is in contradiction to what lRicketj ( Il977ll noted, that there 
is very little flexibility about the Cj va lue in the theo ry. Also, the 
value of C| =5 is in contradiction with lSackeJ l ll974h findings, as 
he obtained Ci = 1.07, incidentally - for the Vela pulsar that we 
also used. Another evidence for C| close to unity was provided by 
ICordes et alj h990h in their analysis of the scintillation and scat¬ 
tering properties of the millisecond pulsar PSR B1937-I-21. Taking 
into account all of the above we suggest that our result should be 
viewed with a dose of scepticism. Nevertheless, similarly to L04 we 
believe that the evidence for larger values of Ci is getting stronger, 
as is the evidence that the value of Cj may be line-of-sight depen¬ 
dant, and that the ISM is in fact highly anisotropic and/or inhomo¬ 
geneous at least along some of the lines-of-sight. 

We based our further analysis on assumption that indeed Ci = 
5, meaning that the middle plot in Fig.|4]is valid. The next step was 
to try to model a relation between at 1 G Hz and the dispersion 
measu re. First we used a relation proposed bv iRamachandran et al.l 
( 1 19971 also used by L04): Tj = A DM^(1 -l- B DM^). Following 
L04 we fixed y = 2.2 (to estimate the deviations from Kolmogorov 
value), which yielded the best fit with: 

rrf(ms) = 2.26 x 10“’ DM’ ’(1 -t 0.00205 DM* ’^). (7) 

The results of this fit are shown in the bottom panel of Fig.|4]as 
a solid blue line. The values of all the parameters are quite different 
from what L04 obtained, although one has to remember that we 
used data for 60 pulsars, while they had only 27 pulsars available, 
and also they assumed Cj = 1. Especially the difference in the value 
of if is striking, as L04 got ( = 23 and earlier lRamachandran et al.l 
d 19971) estimated ( = 2.5. This means that our data indicates a much 


® for the less reliable ratio estimates we got a ratio of 9.7 for B0329+54 
and a ratio of 56 in case of B0823+26, but for both pulsars the values of 
at I GHz from scattering are hugely overestimated due to a very probable 
underestimation of a 


flatter scatter time versus DM dependence, especially in the high- 
DM range. 

Another form o f this dependence we tried was the parabolic 
relation proposed by iBhat et alJ ) l2004ll : adopting this we got the 
best fit corresponding to a relation for scattering at 1 GHz: 

logTrf(ms) = -6.344 -)- 1.467 log DM - 0.509 (log DM)’, (8) 

which is represented by a dashed blue line in the bottom panel of 
Fig-ia As one c an see we got much stronger linear term as for 
[Bhatetal.l ) l2004h it was equal to 0.154 (our value of 1.467 is almost 
10 times g reater), and much lower quadratic term; our 0.509 versus 
1.07 from iBhat et al.l )l2004h . One has to remember however that 
in that paper authors used a different approach to the problem, as 
they used 371 individual scatter time estimates and for the purpose 
of frequency scaling they used the scaling index a as one of the 
fit parameters, and in the fit they assumed it to be the same for 
all the pulsars. Using this method they obtained the average value 
of a = 3.86 ± 0.16, which can be considered remarkably close to 
the average values from our data (see the bottom plot in Fig. IS. 
In our study we used the actual estimates of the frequency scaling 
indices for 60 pulsars (each pulsar is represented by a single mark 
on the plot) to obtain the standardized scatter time value for a given 
object. These indices were different for every pulsar (ranging from 
a = 2.73 to 5.61) and they come from about 300 individual scatter 
time estimates, as each of our a (and thus the 1 GHz tj) values was 
based on at least 3, and up to 9 individual scatter time estimates. 

In the bottom plot if Fig|4]we are also showing the results of 
the previous attempts to find the versus DM relation. The red 
dash-d otted line shows the formula given bv iRamachandran et al.l 
)ll997l) . the dotted line is based on the formula given by L 04, and the 
three-d ot-dashed red line is based upon the formula from iBhat et al.l 
)|2004|) : all these lines were adjusted for the value of C| = 5. With 
this adjustment (which is supposed to work better for the low-DM 
range, where the scintillation based estimates dominate) one can 
see that all the proposed relations basically agree in the low- to 
mid-DM range, especially given the large spread of data. The dis¬ 
crepancies appearing at high-DM values are also probably due to 
the Cl adjustment, but one thing to note is that the formula given 
by L04 (dotted line) seems to be the most similar to our results and 
is showing a similar slope in the high-DM range. This should not 
come as a surprise since our dataset is just an extension of their 
work. 

The differe nces between our approach and the one that 
iBhat et aD | |2004|) used, which in turn influenced the available pul¬ 
sar sample, are probably responsible for the different results of the 
Tj versus DM relation. In both cases however one has to note a large 
scatter present in the data. The source of the scatter may be astro- 
physical in nature (like the strongly line-of-sight dependant nature 
of the scatt ering), as well as m ethod-dependant, and probably both. 
In case of iBhat et al.l )l2004ll one of the sources of an additional 
spread in the data may be the assumption of an uniform value of 
a, what does not agree with observations. In our method the addi¬ 
tional scatter may be present due to using different scaling indices, 
as in such case the method will yield different outcome depending 
on the value of the observing frequency to which we “standardize” 
the amount of scattering. Finally, both methods will suffer due to 
the uncertainties of individual pulse broadening measurements. 

Because of the data scatter we were unable to perform any 
meaningful fit quality testing to check which of the proposed func¬ 
tional forms fits the data better. The values of normalised are 
similar for both models (of the order of 0.44) which does not allow 
to decide which of the fits is more reliable. Large scatter of the data 
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also translates to very low precision of the estimated parameters. 
Also, it is quite possible that the relation has another form, which 
we did not consider. 

Regardless one has to remember that any kind of relation de¬ 
rived from this type of analysis will be only an attempt to aver¬ 
age scattering happening along sometimes vastly different lines-of- 
sight, therefore a large scatter will be always present. The relations 
between the pulse broadening and DM we propose here should be 
at best considered only a very rough estimates. For a given pulsar, 
when calculating the expected scatter time fr om any of these for¬ 
mula s, or the ones propo sed earlier by L04, iRamachandran et^ 
1 19971 or lShat et al.l OOOdI) one has to remember that the actual scat¬ 
tering may be at least one order of magnitude larger or smaller than 
the derived value, depending on the properties of ISM along the 
given line-of-sight. 

Also one has to be aware that this kind of analysis will be 
always prone to errors due to various effects mentioned earlier in 
this paper, such as anisotropies, the inhomogeneity of ISM and our 
lack of knowledge about the distribution of free electrons along the 
line-of-sight. 


5 DISCUSSION AND CONCLUSIONS 

In this paper we performed analysis of the multi-frequency scat¬ 
ter time measurements, following LOl and L04 and focusing on 
pulsars for which the pulse broadening was estimated at multiple 
frequencies. To do so we constructed a database of all the scat¬ 
ter time and decorrelation bandwidth measurements we could find. 
This database was assembled from both our own observations, as 
well as the data available - both in form of pulse profiles as well as 
the data from the literature. It consists of over a 100 pulsars which 
had pulse broadening and/or the decorrelation bandwidth measured 
at multiple frequencies. Over the course of our analysis, and us¬ 
ing similar datasets published previously we were able to gather 
a number of sources for which we were able to reliably ascertain 
the value of the pulse broadening frequency scaling indices which 
are based on scatter time measurements from at least 3 frequencies. 
The largest previous analysis of a similar dataset was done by L04 
(see also references therein) who had the data for 27 objects. Based 
on our data published in Paper 1 we added the (three-frequency 
based) estimates of a for 11 pulsars and some of these were up¬ 
dates of the values published in LOl and L04. As a result of our 
current study we added or updated the a and (at the standard fre¬ 
quency of 1 GHz) values for 33 sources, which yielded a total of 60 
pulsars with reliable estimates of the scattering frequency scaling 
indices, over two times the number L04 had. 

A large portion of resulting values of the a indices does not 
conform with the predictions based on the most commonly used 
models of the scattering phenomenon. The simplest model of a 
scattering with a thin screen geometry and Kolmogorov-type den¬ 
sity fluctuation within the screen predicts the pulse broadening 
spectral slope of 4.4. Using other fluctuation models one can get 
expected values of a down to 4.0, and applying different geome¬ 
tries of ISM along the line-of-sight one can predict a values lower 
than 4.0. 

We extensively discussed the phenomena that influence 
the measured value of a in our previous paper (Paper 1, 
iLewandowski et aT]|2013h . and virtually all the arguments pointed 
there still hold. There are a few reasons for which the measured 
pulse broadening spectral slope may be lower. The main is probably 
geometry of ISM, including the possibility of scattering on multiple 


screens (see also LOl and L04). If the actual spe ctrum of the fluc¬ 
tuation density different from Kolmogorov’s (see lGuptil2000l for a 
brief summary) it will also play its role in the general picture. An- 
other source of disc repancy may be the inner scale effects (see also 
iRickett et al.ll2001^ which can affect and flatten the spectral slope 
below certain observing frequency. Finally, the anisotropy and the 
inhomogeneity of the turbulence may also influence the observed 
scatter time spectra, both directly, as well as indirectly - both of 
these effects may cause the scattering parameters to vary in time, 
hence they may make it impossible to create a proper spectra of 
pulse broadening. 

Our current analysis mainly due to its completeness (with re¬ 
gard to the previous similar attempts) allows us also to study the de¬ 
pendence of basic scattering parameters with dispersion measure. 
As we mentioned above, prior t o our work the most comple te anal¬ 
ysis of this kind was made by iLdhmer et alj ( l200lLl2003l . Based 
on the data they had available they postulated, that the scatter time 
spectral index a stays with good agreement with the theoretical pre¬ 
dictions for pulsars with DM< 300 pc cm“^, while for objects with 
larger DM the slope is much lower (a = 3.48 on average). 

In Paper I we modified this conclusion only slightly, but with 
the new or updated estimates of a we obtained (especially for pul¬ 
sars with low- and medium-DMs) it seems that there is no clear 
division between the low-DM and high-DM pulsars. The devia¬ 
tions from the theoretically allowed range of scatter time scaling 
indices (4.0 < a < 4.4) start from the lowest DM values for se¬ 
lected sources. If anything, the range of the deviations seem to in¬ 
crease with higher dispersion measures. Also, one has to note that 
not all mid- and high-DM objects deviate from theoretical predic¬ 
tions, and it seems that it is only the probability of deterring from 
theory that increases for larger DMs. 

Collecting a estimates for 60 pulsars allowed us to average the 
values over various DM intervals (see the bottom plot of Fig.|2]l. 
The result was that in a broad range of dispersion measures the 
average value of the spectral scaling indices stays almost the same 
and only slightly below a = 4.0, and that is despite the fact that for 
individual sources the value of the scaling index can vary a lot. Only 
for DM>500 pc cm“^ we noted a slight deviation of the average 
a value, but becomes statistically significant for very large DMs 
(>700 pc cm^'^). The average value of the scaling index from all 
of our data is a = 3.89, for pulsars with DM<500 pc cm“^ the 
average is a = 3.95 and for sources with DM>500 pc cm“^ the 
average value is a = 3.49 (based on 8 pulsars). 

As it was previously known the simplistic thin-screen model 
can not explain the scattering properties for very distant objects. Es¬ 
pecially with the possibility of multiple scatterings along the line- 
of-sight one could expect, that the number of such occurrences with 
increase for more distant pulsars, however there is no guarantee for 
that, as most of the astrophysical objects that could play the role 
of a thin scattering screen (such as dense HII regions) can be rel¬ 
atively small. So, while the chance that the line of sight will cut 
such screen will increase with the distance, there may still be some 
selected directions in which that will not happen. The probability 
of undisturbed line-of-sight however should rapidly decrease with 
distance, and seemingly this is what we observe. 

We also analysed the actual amount of scattering visible in 
pulsars from our sam ple versus the disp ersion measure. Similarly 
to L04, and contrary to iBhat et al.l ( I 2 OO 4 I 1 for each individual source 
we used the actual scatter time frequency scaling law to estimate 
the amount of pulse broadening at the frequency of 1 GHz. The re¬ 
sults are shown in Fig. |4] To bind the scatter time estimates from 
the actual pulse broadening with the values obtained from decor- 
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relation bandwidth measurements we decided to adjust the value 
of the Cl constant which appears in the relation connecting these 
two quantities (see Eq.[3j- Based on measurements of of the Vela 
pulsar, and B1933+16 we used the value of Ci = 5, although one 
may argue that the evidence for such value is not very strong. To 
describe the Tj(1 GHz) versus DM dep endence we tried to model i t 
using the formulae proposed earlier by|Ramachandran_e£^ l ll997ll 
as well as the parabolic relation used by Bhat et alj l l2004h . In both 
cases the resulting formulas are seemingly quite different from the 
ones published previously (when one considers the fitted parame¬ 
ters), however upon inspecting the plots they seem in agreement 
with each other, at least in the low- to mid-DM range. If anything, 
our Tj versus DM relations look more linear (on the log-log plot) 
than in the previous attempts to recreate it. This also means that for 
the mid-DM pulsars we (on average) expect more scattering, while 
for high-DM objects the rise of the scattering strength will not be 
as sharp as it was previously predicted. 

Our database of multi-frequency scatter time measurements 
and the frequency scaling index estimates presents a significant ex¬ 
tension of the previous similar attempts, there is however still a lot 
of possibility for improvement. The frequency coverage can be, and 
probably will be better in the nearest future, given the fact that the 
advances in technology will allow us to use wide-band receiver sys¬ 
tems even at radio frequencies below 1 GHz (one of such examples 
is the GMRT upgrade which is currently happening). It will defi¬ 
nitely be improved even more for the low-DM pulsars with the help 
of the low-frequency telescopes such as LOFAR. The future looks 
bright, which is really good, because - as our analysis presented 
here shows - we still do not know everything about the scattering 
of pulsar radiation in our Galaxy. 

Finally, want to emphasize probably the most important out¬ 
come of our analysis: it appears that the properties of scattering of 
radio waves in the Milky Way are much more line-of-sight depen¬ 
dant than it was previously thought, even for more distant objects. 
Thus any attempts to generalize the behaviour of scattering may 
prove to be futile, and large scatter observed in the scattering fre¬ 
quency scaling index, or the scattering strength, versus the disper¬ 
sion measure plots will be always present. 
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APPENDIX: MULTI FREQUENCY SCATTER TIME 
MEASUREMENTS DATABASE 

In this appendix we present the results of the scatter time mea¬ 
surements, both our own as well as the ones found in the literature, 
that were included our multi-frequency scatter time database, along 
with the data published earlier in Lewandowski et al. (2013). We 
present the data for all the pulsars from our studies for which we 
obtained new measurements, either from our own observations or 
from the archived profiles. 


Below is the list of references to the table: 

(1) Alukar, Bobra & Slee (1986) 

(2) our analysis based on the profiles from the EPN Database 

(3) Rickett (1977) 

(4) Manchester, D’Amico & Tuohy (1985) 

(5) Rankin (1981) 

(6) Cordes, Weisberg & Boriakoff (1985) 

(7) Gullahorn & Rankin (1978) 

(8) Slee, Dulk & Otrupcek (1980) 

(9) Lewandowski et al. (2013) 

(10) Ldhmer et al. (2001) 

(11) Abies, Komesaroff & Hamilton (1970) 

(12) our analysis based on the pulsar profiles from the ATNF 
Parkes Database 

(13) Kuzmin & Losovsky (2007) 

(14) our analysis of new observations performed at GMRT 

(15) Staelin & Shutton (1970) 

(16) Ulyanov, Seredkina & Shevtsova (2012) 


Table 3. Full listing of the scatter time measurements 


Pulsar 

/obslMHz] 

Tj [ms] 

Reference* 

B0136+57 

102 

25 6 

13 


408 

1.43±0.1 

2 


610 

1.23±0.04 

2 


925 

1.5±0.08 

2 


1408 

1.55±0.05 

2 


1642 

1.56±0.88 

2 

B0226+70 

111 

13±2 

13 


1408 

0.3 

2 


1420 

1.21±0.57 

2 

B0320+39 

40 

130±40 

13 


102 

0.4 

2 


408 

1.02 

2 



1.44±0.14 

2 



0.39 

2 



0.5 

2 


610 

0.47 

2 



0.95 

2 



0.43±0.04 

2 



4.39 

2 



4.92 

2 


910 

15.72±3.51 

2 


1400 

2.68±0.42 

2 


1408 

20.48±1.52 

2 


4850 

22.28±7.55 

2 

B0339+53 

111 

25±5 

13 


606 

1.45±0.64 

2 



1.54 

2 


1400 

0.0057±0.001 

2 

B0353-t52 

111 

45±10 

13 


1420 

3.32±0.05 

2 

B0402+61 

111 

28±6 

13 


400 

12.36±1.01 

2 


800 

0.19±0.11 

2 


925 

0.16±0.07 

2 


1330 

0.31 ±0.0002 

2 


1408 

0.11 

2 

B0450-18 

111 

5.5±3 

13 


150 

5.8±1,6 

1 


408 

0.57 

2 


610 

0.52±0.22 

2 


925 

0.473±0.01 

2 


1404 

0.27±0.06 

2 


1642 

0.29 

2 


(continued on 

following pages) 



* full list of references at the beginning of the Appendix 
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Table|3l Scatter time measurements (continued). Table|3l Scatter time measurements (continued). 


Pulsar 

/ob.[MHz] 

Tj [ms] 

Reference* 

Pulsar 

/obs[MHz] 

Trf [ms] 

Reference* 

B0531+21 

44 

300±80 

13 

B0809+74 

23 

3 

16 


63 

100±20 

13 



10 

16 


111 

10±3 

13 


40 

4 

16 


318 

0.23 

15 


44 

3.9+1 

13 


2695 

0.002 

2 


62 

0.5 

16 






63 

0.6 0.2 

13 

J0533+0402 

102 

70±30 

13 


102 

0.2 

16 


370 

6.6+1.5 

2 


410 

12.83+0.44 

2 






606 

8.37+0.67 

2 

B0540+32 

111 

15±5 

13 


925 

0.05+0.01 

2 


610 

3.8±0.08 

2 


1408 

0.63+0.02 

2 

B0559-05 

111 

28±10 

13 

B0823+26 

44 

12+3 

13 


400 

6.68+0.8 

2 


63 

4+2 

13 


800 

6.98+0.56 

2 


400 

20 

3 






408 

0.23+0.02 

2 

B0611+22 

111 

40+10 

13 


1055 

0.009+0.001 

2 


925 

2.37+0.23 

2 






1408 

0.98+0.2 

2 

B0833-45 

160 

621+109 

1 


1420 

0.21+0.09 

2 



170+30 

8 


4850 

0.23 

2 


250 

20 

11 

B0621-04 

111 

45+4 

13 


300 

9.4 

3 


610 

1.95+0.92 

2 



12 

11 



1.99 0.93 

2 


318 

6 

2 


1400 

1.74+0.27 

2 


408 

4.7 

11 



1.76+0.27 

2 


636 

2 

11 


1408 

1.71+0.59 

2 


1392.5 

2.25+0.08 

2 

B0626+24 

111 

55+30 

13 


1420 

1.12 

11 


1408 

0.99+0.2 

2 

B0835-41 

160 

24+5.5 

1 







30+6 

8 

B0740-28 

102 

22+5.3 

13 


685 

0.86+0.08 

12 


160 

24.5+2.8 

1 


1590 

0.83+0.87 

12 



21+3 

8 






408 

5.34+0.39 

2 

B0839-53 

410 

20+2.4 

1 


422 

1.09+0.2 

2 


660 

2.14+0.24 

2 


660 

1.41+0.1 

12 


1440 

0.24+3.8 

2 


1390 

1.16+0.17 

12 






1518 

0.05+0.01 

12 

B0919+06 

44 

25+15 

13 






63 

9+5 

13 

B0808-47 

410 

33+1.5 

1 


80 

15+2 

8 


660 

7.18+0.33 

2 


111 

1.5+0.8 

13 


1326.75 

1.94+0.7 

12 


408 

0.09+0.01 

2 


1382.5 

1.46 

2 


610 

0.32 

2 


1440 

5.14+0.94 

2 


925 

0.04 

2 






1420 

0.44 

2 

* full list of references at the beginning of the Appendix 


1710 

0.26 

2 





B1054-62 

410 

23+2.5 

1 






1382.5 

8.52+0.35 

2 






1560 

8.46+0.24 

2 


* full list of references at the beginning of the Appendix 
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Multi-frequency pulsar scatter time estimates 


Table|3l Scatter time measurements (continued). Table|3l Scatter time measurements (continued). 


Pulsar 

/obs[MHz] 

Trf [ms] 

Reference* 

Pulsar 

/ob.[MHz] 

Trf [ms] 

Reference* 

B1114-41 

160 

7.4±2.9 

1 

B1558-50 

410 

20.8±2.2 

1 


658 

1.36±0.62 

2 


1374 

0.55 

12 


660 

0.89±0.09 

2 


1518 

1.88±0.78 

12 


1318.5 

0.08±0.01 

12 


1693 

0.1 

12 

B1154-62 

410 

27±1.7 

1 

B1620-42 

410 

36±4.7 

1 


685 

3.86±0.26 

12 


1440 

3.17±1.93 

2 


1327 

0.12±0.01 

12 






1590 

0.07±0.01 

12 

B1641-45 

610 

73.4±4.8 

9 






750 

40 

3 

B1302-64 

410 

25±2.1 

1 


1060 

8.8±0.4 

9 


1382.5 

0.58±0.06 

2 


4850 

2.5±0.4 

9 

B1323-58 

410 

77±15 

1 

B1642-03 

63 

7.2±1.2 

13 


750 

30±0.27 

3 


80 

4.4±0.11 

8 


1431 

2.21±0.23 

12 


111 

1.5±0.5 

13 


1560 

0.08±0.02 

2 


408 

0.9±0.02 

2 


1676 

1.38±0.38 

12 









J1705-3950 

610 

35±9 

14 

B1323-62 

410 

185±35 

1 


1374 

3.89±0.23 

12 


685 

19.42±0.15 

12 






750 

30±0.27 

3 

J1723-3659 

610 

17.2±2.5 

14 


1318 

0.4±0.04 

12 


1374 

1.3±0.32 

12 


1347 

0.15 

12 


1518 

1.03±0.63 

12 


1382.5 

0.5 

2 






1490.2 

0.58 

2 

B1737-39 

410 

11 ±0.8 

1 






658 

4±0.09 

2 

B1356-60 

410 

24±1.7 

1 


660 

4.39±0.13 

12 


660 

5.05±0.05 

12 






1374 

0.23 

12 

J1739-3023 

610 

0.28 

14 


1693 

0.57±0.04 

12 


1374 

0.03±1.37 

12 


2350 

0.01 

12 


1518 

0.04±4.15 

12 

B15084-55 

44 

15±4 

13 

J1744-3130 

610 

11.8±6.0 

14 


63 

6.4±4 

13 






408 

0.62±0.24 

2 

B1745-12 

111 

55±24 

13 


610 

0.18±0.002 

2 


610 

5.06±0.2 

2 



0.14 

2 


1382.5 

7.23±0.56 

2 



0.35 

2 






925 

0.5 

2 

B1749-28 

80 

24.9±2 

1 


1408 

0.63±0.002 

2 



40±5 

8 


1642 

0.35 

2 


160 

2.9±0.1 

1 


4750 

5.26±1.9 

2 



3.1±0.5 

8 






408 

1.58±0.01 

2 

B1557-50 

325 

411±150 

9 


610 

1.26±0.05 

2 


610 

43±3 

9 


925 

1.83±0.058 

2 


750 

16±0.25 

3 


1408 

2.02±0.035 

2 


1060 

5±0.6 

9 


1642 

0.31±0.09 

2 


1400 

1.98-1-0.35 

9 






1600 

0.97±0.24 

9 

* full list of references at the beginning of the Appendix 


* full list of references at the beginning of the Appendix 
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Table|3l Scatter time measurements (continued). Tables Scatter time measurements (continued). 


Pulsar 

/ob.[MHz] 

Td [ms] 

Reference* 

Pulsar 

/ob.[MHz] 

Td [ms] 

Reference* 

B1758-03 

111 

55±30 

13 

B1846-06 

111 

100+40 

13 


610 

0.12+1.8 

2 


610 

3.42+0.5 

2 

B1758-23 

1274.75 

130.53+5.38 

12 



3.42+0.45 

2 


1374 

102.49+1.13 

12 



3.39+0.43 

2 


1400 

111+19 

10 


657 

3.47+0.83 

2 



99+19 

4 






1421 

83.22+3.92 

12 

J1852-0635 

610 

15.1+3.6 

14 


1518 

74.32+1.27 

12 


1374 

0.95+0.0004 

12 


1642 

55+10 

10 







51 + 10 

9 

B1853+01 

111 

26+10 

13 


2263 

17.87+0.84 

12 


606 

0.97+0.44 

2 


2600 

0.75+0.34 

9 






4850 

0.23+0.08 

9 

B1859+03 

320 

187 

3 






325 

241+51 

10 

B1802+03 

102 

33+10 

13 


400 

83.93+2.07 

2 






408 

89+9 

10 

B1818-04 

102 

60+40 

13 


410 

68.3+2.29 

2 


160 

50.3+5 

1 


430 

65+11 

5 



65+13 

8 


606 

12.98+0.18 

2 


410 

4.74+0.07 

2 


610 

12.6+0.9 

10 


610 

2.82+0.03 

2 


925 

1.92+0.47 

2 


925 

3.47+0.12 

2 


1408 

0.13+0.01 

2 


1408 

3.28+0.09 

2 


1420 

2+0.13 

2 


1642 

3.47+0.16 

2 


1642 

0.06+0.01 

2 

B1821-19 

325 

57+23 

10 

B1900+01 

243 

176+58 

10 


408 

22+4 

10 


325 

34+2 

10 


410 

17+1.3 

1 


408 

14.7+1.3 

10 


610 

4.6+0.2 

10 


410 

13.04+0.39 

2 


925 

0.92+0.04 

2 


430 

11+2 

5 


1642 

0.00057+0.00018 

2 


610 

3.37+0.07 

2 

B1831-04 

102 

50+25 

13 


925 

1.02+0.12 

2 






1408 

0.71+0.07 

2 

B1839-04 

111 

150+15 

13 









B1907+02 

111 

140+50 

13 

J1841-0345 

610 

15.3+5.2 

14 


160 

27.9+16 

1 






320 

3+1 

6 

B1844-04 

160 

63.5+30 

1 


410 

2.84+0.79 

2 


408 

15.4+0.53 

2 


610 

3.66+0.34 

2 


610 

0.77+0.08 

2 


1408 

5.76+0.68 

2 


925 

0.34+0.04 

2 






1408 

0.04+0.01 

2 

B1907+10 

160 

26.5+8.1 

1 







27+7 

8 

B1845-01 

235 

262+181 

10 


320 

2.2+0.7 

6 


408 

54+8 

10 


410 

0.13+0.05 

2 


410 

42.6+3.4 

2 






430 

42+6 

5 

B1911-04 

102 

35+15 

13 


610 

16+0.19 

2 


160 

16.7+1.8 

1 



17+1.4 

10 



32+5 

8 


925 

5.71+0.36 

2 


408 

0.861+0.6 

2 


1408 

5.95+0.32 

2 


610 

0.72+0.81 

2 


1642 

8.41+0.53 

2 


925 

0.47+0.96 

2 


4750 

0.16+0.05 

2 


1408 

0.12+0.004 

2 






1642 

0.26+1.96 

2 

* full list of references at the beginning of the Appendix 
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0.09 

2 


* full list of references at the beginning of the Appendix 
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Multi-frequency pulsar scatter time estimates 


Table|3l Scatter time measurements (continued). 


Tables Scatter time measurements (continued). 


Pulsar 

/obJMHz] 

Trf [ms] 

Reference* 

B1913+10 

408 

23.87+0.77 

2 


430 

14.6 

7 


610 

5.41+0.15 

2 


925 

1.12+0.37 

2 


1408 

1.56+0.011 

2 


1420 

2.61+0.29 

2 


1642 

0.086+0.052 

2 

B1914+13 

320 

5.5+2.5 

6 


610 

1.18+0.27 

2 


1408 

0.39+0.01 

2 

B1915+13 

102 

40+20 

13 


160 

11.7+1.9 

1 



12+3 

8 


408 

1.66+0.07 

2 


925 

1.76+0.07 

2 


1642 

1.22+0.15 

2 

B1919+20 

430 

11.7 

5 


606 

0.13 

2 

B1919+21 

44 

14.9+3.6 

13 


63 

4.6+3.2 

13 


111 

0.6+0.4 

13 


408 

0.33+0.002 

2 


610 

0.35+0.001 

2 

B1920+21 

102.75 

26 

2 


160 

96.8+50 

1 


243 

4.4+1.5 

10 


410 

5.11+0.95 

2 


925 

8.63+1.37 

2 

B1929+20 

410 

10.58+0.54 

2 


430 

4.83 

7 


925 

0.007 

2 


1408 

0.04+0.04 

2 

B1933+15 

430 

5.8 

7 

B1933+16 

110 

67 

3 


111 

50+15 

13 


160 

21.7+1.6 

1 



25+4 

8 


243 

4.6+0.2 

10 


325 

1.8+0.1 

10 


410 

1.62+0.22 

2 


610 

0.51+0.33 

2 


925 

0.05+0.01 

2 


1408 

0.02+0.01 

2 


* full list of references at the beginning of the Appendix 


Pulsar 

/obs[MHz] 

Trf [ms] 

Reference* 

B1940-12 

160 

4.4+1.1 

1 


410 

0.93 

2 


610 

1.05+0.07 

2 


1318 

1.66+3.58 

12 


1642 

2.26+2.69 

2 

B1943+18 

430 

15+5 

6 

B1946+35 

111 

45+10 

13 


320 

32 

3 


408 

16+0.12 

2 


430 

15+3 

5 


610 

3.71+0.1 

2 


925 

0.6+0.57 

2 


1642 

0.4+0.06 

2 

B1953+50 

111 

45+10 

13 


408 

1.55+0.25 

2 


925 

0.18+0.09 

2 

B2002+31 

243 

27+13 

10 


320 

10 

3 


325 

5.3+1.2 

10 


408 

2+1.4 

10 



5.51+0.76 

2 


610 

2.62+1.09 

2 


925 

0.033+0.049 

2 


1408 

5.33+0.41 

2 


1642 

4.13+1.08 

2 

B2020+28 

63 

2.3+1 

13 


111 

0.5+0.3 

13 


410 

0.83+0.08 

2 

B2027+37 

111 

132+10 

13 


610 

3.94+0.86 

2 



0.5 

2 



3.85+0.81 

2 


1408 

7.15+0.72 

2 



7.16+0.74 

2 


1420 

6.02+0.48 

2 

B2053+36 

111 

74+25 

13 


320 

7.7+1.5 

6 


430 

2.2+0.4 

6 


925 

0.13+0.01 

2 

B2111+46 

111 

120+40 

13 


408 

1.5+2.71 

2 


610 

5.19+0.71 

2 



2.34+0.23 

2 


1408 

2.67+0.39 

2 


* full list of references at the beginning of the Appendix 
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Table|3l Scatter time measurements (continued). 


Pulsar 

/obJMHz] 

Trf [ms] 

Reference* 

B2113+14 

102 

25±10 

13 


1408 

0.11 ±0.03 

2 

B2217+47 

44 

70±20 

13 


63 

12±9 

13 


111 

3.2±0.3 

13 


410 

0.42±0.23 

2 



0.37±1.94 

2 


610 

0.47±4.69 

2 


1408 

0.5 

2 

B2303+46 

325 

1.06±0.0005 

14 


610 

0.65±0.03 

14 


1060 

0.35 

14 

B2303+30 

44 

300±100 

13 


63 

110±20 

13 


111 

13±3 

13 


160 

9.9±3.6 

1 


410 

1.47 

2 


610 

1.66 

2 


1408 

0.631 

2 

B2351+61 

111 

35±5 

13 


* full list of references at the beginning of the Appendix 
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